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Abstract 


2 This paper gives an overview of August 2004 through July 2009 upper tropospheric 

3 (UT) water vapor (H2O) and ice water content (IWC) from the Aura Microwave Limb 

4 Sounder (MLS) and comparisons with outputs from the NASA Goddard Earth Observing 

5 System Version 5 (GEOS-5) data assimilation system. Both MLS and GEOS-5 show that 

6 high values of FLO and IWC at 215 to 147 hPa are associated with areas of deep 

7 convection. They exhibit good (within -15%) agreement in IWC at these altitudes, but 

8 GEOS-5 HiO is -50% (215 hPa) to -30% (147 hPa) larger than MLS, possibility due to 

9 its higher temperatures at these altitudes. GOES-5 produces a weaker intertropical 

10 convergence zone than MLS, while a seasonally-migrating band of tropical deep 

11 convection is clearly evident in both the MLS and GEOS-5 UT FLO and IWC. MLS and 

12 GEOS-5 both show spatial anti-correlation between IWC and H 2 O at 100 hPa, where less 

13 FLO is associated with low temperatures in regions of tropical convection. At 100 hPa, 

14 GEOS-5 produces 50% less IWC and 15% less FLO in the tropics, and -20% more H 2 O 

15 in the extra-tropics, than does MLS. Behavior of the 100 hPa FLO, which exhibits a 

16 quasi-biennial oscillation, appears consistent with it being controlled by temperature. The 

17 seasonal cycle in the vertical transport of tropical mean FLO from -147 hPa to -10 hPa 

18 appears much stronger in MLS than in GEOS-5. The UT IWC and FLO interannual 

19 variations, from both MLS and GEOS-5, show clear imprints of the El Nino-Southern 

20 Oscillation. 
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1. Introduction 


2 Upper-tropospheric (UT) water vapor (H 2 O) and clouds play important roles in 

3 regulating Earth’s climate, producing feedbacks to climate forcings by increasing 

4 greenhouse gases. HoO is the primary natural atmospheric greenhouse gas, trapping some 

5 of the outgoing longwave radiation (OLR) that would otherwise be emitted to space. The 

6 increase of UT fUO with sea surface temperature (SST) provides a strong positive 

7 feedback in response to surface temperature increases that can be caused by increasing 

8 anthropogenic greenhouse gases [e.g. Held and Soden, 2000; Su el al. 2006a], Udelhofen 

9 and Hartmann [1995] showed that OLR is mostly sensitive to UT relative humidity 

10 changes above 400 hPa. Climate models indicate that UT H 2 O could increase -200% by 

11 the end of the 21 st century, compared to a -20% increase in lower tropospheric H 2 O 

12 [Soden et al., 2005]. This UT amplification underscores the importance of monitoring 

13 and quantifying UT H 2 O variability. 

14 Clouds in the UT tend to have a net warming effect, as their cold tops result in low 

15 OLR [ Stephens , 1990; Su et al, 2008]. The occurrence of UT clouds is closely related to 

16 UT humidity [ Udelhofen and Hartmann, 1995; Soden and Fu, 1995; Su et al., 2006a], 

17 The variation of UT cloud amount with SST and the resulting potential climate feedback 

18 have been a subject of debate [ Lindzen et al., 2001; Lin et al, 2002; Hartmann and 

19 Michel sen, 2002; Su et al, 2008], The UT cloud radiative heating also influences 

20 transport from the troposphere to the stratosphere [e.g. Cord et al, 2006; Hartmann et 

21 al., 2001]. 

22 The Microwave Limb Sounder (MLS) on board the Aura satellite, launched on July 

23 15, 2004, provides simultaneous global measurements of UT H 2 O, cloud ice water 

24 content (IWC), temperature (T), and several trace gases [Waters et al. 2006]. Li et al. 


2 



1 [2005; 2007] compared Aura MLS IWC measurements with European Centre for 

2 Medium-range Weather Forecast (ECMWF) analyses and forecasts, and with other state - 

3 of-the-art climate model simulations, and found differences as large as factor of 10 

4 between models and observations. These helped promote modifications to the ECMWF 

5 model cloud microphysics that resulted in significant improvement [Waliser, et al., 

6 2009]. Su et al. [2006] found differences between models and observations of up to a 

7 factor of 4 in the relationships among UT PLO, IWC, and SST. Read et al. [2008], using 

8 MLS PLO and CO measurements, estimated the relative influences of convection, 

9 “freeze-drying” and extra-tropical mixing on the amount of H 2 O entering the 

10 stratosphere. 

11 This paper presents an overview of the global distributions and temporal variations 

12 for UT IWC and PLO as seen by MLS from August 2004 through July 2009 (the five- 

13 year period for which data are currently available), spatial correlations with deep 

14 convection and temperature, and comparisons with output from the Goddard Earth 

15 Observing System Version 5 (GEOS-5) data assimilation system. Section 2 describes the 

16 datasets, section 3 presents spatial distributions, and section 4 presents temporal 

17 variations. Section 5 focuses on the UT response to the El Nino-Southern Oscillation 

18 (ENSO). Section 6 gives conclusions and discussion. 

19 2. Data 

20 2.1. Aura MLS UT Water Vapor, Cloud, and Temperature Measurements 

21 We use MLS Version 2.2 (V2.2) Level 2 [Livesey et al., 2007] H 2 O, IWC and T 

22 datasets, whose validations are described by Read et al. [2007], Wu et al. [2008], and 

23 Schwartz et al. [2008], respectively. MLS measures -3500 vertical profiles per day along 

24 a sun-synchronous suborbital track having equatorial crossings at 1:40 PM and 1:40 AM 
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1 local solar times. The Level 2 data are produced on pressure surfaces (12 surfaces per 

2 decade) from 316 to 0.1 hPa, with IWC having a limited useful range of 215 to 83 hPa. 

3 These data have a vertical resolution of ~3-4 km, and horizontal resolutions of ~7 km 

4 across-track and -200-300 km along-track. Estimated measurement accuracies are 20% 

5 for FLO. 2 K for temperature (which, in cloudy regions, has a known low bias of -2 K at 

6 215 hPa and -0.5 K at 147 hPa [ Schwartz et al. 2008]), and a factor of two for IWC. 

7 MLS measurements are generally not degraded by the presence of clouds and 

8 aerosols, whose particle sizes are typically much smaller than the measurement 

9 wavelengths. Very thick clouds (IWC > -50 mg/m3) can degrade the temperature and 

10 some species measurements [Wu et al. 2008], but the retrieval algorithms [Livesey et al., 

11 2006] flag such measurements and they are not used here. See Wu and Jiang [2004] for 

12 details on the identification and quantification of cloud- affected radiances and the IWC 

13 retrieval. 

14 2.2. GEOS-5 Meteorological Products 

15 Meteorological datasets produced by the GEOS Versions 5.1.0 and 5.2.0 data 

16 assimilation systems are used in this study. Rienecker et al. [2008] describe the 

17 meteorological analysis, which uses a three-dimensional variational (3D-Var) approach 

18 [Sasaki, 1970]. GEOS-5. 1.0 ran in near-real time between 17 October 2007 and 14 

19 August 2008; it was also used to retroactively analyze the period from October 2003, 

20 before the Aura launch, until October 2007. GEOS-5. 1.0 was replaced by GEOS-5. 2.0 on 

21 14 August 2008. Both versions, collectively referred to as GEOS-5, produce analyses, 

22 forecasts and assimilated fields on a 72-layer grid, extending from the surface to 0.01 

23 hPa, with a 0.5° x 0.67° latitude-longitude resolution. Vertical resolution is -1.5 km in 
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1 the UT. Differences between the two versions of GEOS-5 will be mentioned as necessary 

2 in the presentation of results. 

3 The GEOS-5 analyses are “snapshots” of the atmospheric state produced four times 

4 daily (at 00Z, 06Z, 12Z, and 18Z) using optimal combinations of model forecasts and 

5 many observations [Rienecker et al., 2008] via the Grid-point Statistical Interpolation 

6 (GSI) technique of Wu et al. [2002]. The assimilated fields are a continuous time series 

7 produced using the GEOS-5 atmospheric general circulation model (AGCM), in which an 

8 additional forcing term is added to the momentum, thermodynamic, moisture and ozone 

9 equations, following Bloom et al. [1996]. This “incremental analysis update” (IAU) 

10 forcing is computed from the difference between analysis and AGCM forecast at the 

11 analysis times, then added as a forcing tendency that remains constant in six-hour 

12 segments that straddle the analysis times. The assimilated data are these AGCM fields 

13 that are constrained by the analyses and contain all information derived from the model, 

14 such as cloud and radiation fields, in addition to the analyzed variables [Rienecker et al., 

15 2008], 

1 6 The importance of transport to the moisture budget and the fact that all clouds in the 

17 GEOS-5 assimilations depend strongly on the AGCM require that some details of the 

18 model be mentioned for understanding the products. The GEOS-5 AGCM is coded 

19 flexibly, but used in particular configurations (spatial resolution and physical parameter 

20 settings) in each version of the assimilation system. Adiabatic transport is computed 

21 using the “finite- volume dynamical core” [Lin, 2004] with a quasi-Lagrangian vertical 

22 coordinate, followed by remapping to the standard 72-layer hybrid grid on which 

23 physical tendencies are computed every 30 minutes. The model includes prognostic 

24 equations for large-scale gases, liquid (condensate) and ice (anvil-type) water, with 
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1 consideration of sub-grid convective contributions to the large-scale liquid and ice 

2 phases. Convection is computed using an adaptation of the Relaxed Arakawa-Schubert 

3 (RAS) convection code [Moorthi and Suarez, 1992], with modifications based on Sud 

4 and Walker [1999] as described in Bacmeister et al. [2006]. RAS considers a sequence of 

5 detraining convective plumes extending between cloud base (set as a fixed layer in 

6 GEOS-5, but inherently adaptable in RAS) and each layer below the tropopause region 

7 (close to lOOhPa); each plume produces detraining mass and cloud condensate at each 

8 layer and also modifies the environmental meteorological (temperature, moisture, wind) 

9 profiles felt by the next plume. The large-scale cloud condensate scheme, based on 

10 probability distribution functions (PDFs) of the moisture field assumed inside a grid box, 

11 incorporates changes to condensate and anvil clouds obtained from RAS, then computes 

12 new sources for the anvil cloud (freezing of existing condensate) and new partitioning of 

13 condensate, before computing loss due to evaporation, auto-conversion of liquid or 

14 mixed-phase condensate, sedimentation of frozen condensate, and accretion of 

15 condensate by falling precipitation. Details of these processes are given in Rienecker et 

16 al. [2008], 

17 Atmospheric moisture in GEOS-5 is analyzed in the form of relative humidity along 

18 with other analysis variables including the stream function, the unbalanced part of 

19 velocity potential, temperature, surface pressure, ozone, cloud liquid and ice water, and 

20 regression coefficients for radiance bias correction. The optimal analysis is obtained by 

21 finding the best fit to the six-hour forecast field and observations while minimizing the 

22 cost function. Various observation types such as radiosondes and radiances from the 

23 Microwave Humidity Sounder (MHS), the Special Sensor Microwave/Imager (SSM/I) 

24 and the Atmospheric Infrared Sounder (AIRS) provide information to constrain the 
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1 moisture fields. AIRS, in particular, gives information on vertical structure for 

2 atmospheric temperature and moisture due to its sounding capability. A total of 152 

3 spectral channels from AIRS are currently assimilated in GEOS-5; these are selected 

4 from the 281-channel “NWP subset” of AIRS radiance measurements. While most of the 

5 AIRS channels are subject to water vapor absorption, it is most significant in the infrared 

6 portion of the spectrum from 6.20pm to 8.22pm. At present, 49 of the water vapor 

7 channels from this band are being assimilated. These water vapor absorption channels 

8 peak at different pressure levels in the troposphere, providing information on the vertical 

9 distribution of moisture for the analysis. The observation error covariance matrix for 

10 radiance data is assumed to be diagonal, i.e., possible inter-channel correlations are 

11 neglected. The error values assigned to the water vapor channels are larger than those 

12 assigned to temperature channels in order to account for the possibility of inter-channel 

13 error correlations, the effects of undetected residual cloud, and the non-linear nature of 

14 the moisture channels that is not accounted for in the formulation of Jacobians in the 3D- 

15 Var analysis. 

16 For comparison with MLS, the GEOS-5 data are interpolated onto the MLS 

17 measurement locations in both space and time. Previous studies [e.g., Li et cil. 2007; Su et 

18 al. 2006b] have shown that such interpolation is particularly important because of 

19 potential artifacts that arise from incomplete sampling of the diurnal cycle by polar- 

20 orbiting satellites. For horizontal sampling, GEOS-5 data are collocated with MLS data 

21 by averaging the data in boxes of 3° along the track and 1° across the track centered on 

22 the MLS measurement locations (approximately matching the MLS footprints). 

23 Vertically, the respective MLS averaging kernels [Read et al. 2007, Schwartz et al. 2008] 

24 are applied to the GEOS-5 H 2 O mixing ratio and temperature products. The GEOS-5 
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1 IWC data are averaged in vertical boxes of ~3.5 km centered on MLS data points to 

2 mimic the MLS IWC vertical resolution [Wu et al. 2008]. 

3 3. Global distribution of UT H 2 0 and IWC 

4 3.1. Annual mean maps 

5 Figure 1 shows four-year mean (January 2005 to December 2008) annual IWC, FLO 

6 and T maps at three pressure levels (100, 147, and 215 hPa), from both MLS observations 

_2 

7 and GEOS-5 analyses. Contours enclosing GEOS-5 OLR of 240 Wm “ or less (indicating 

8 regions of deep convection), and potential vorticity contours of 3.5 x 10 _6 Km 2 kg _1 sec _1 

9 (PV3.5, indicating the poleward edge of the dynamical tropopause), are superimposed. 

10 Both MLS and GEOS-5 data show that at 215 and 147 hPa, large IWC and FLO and low 

11 OLR are collocated in the tropical western Pacific, west central Africa and northern 

12 South America. The PV3.5 contour generally encloses the large IWC and FLO values, 

13 supporting the notion [e.g. Wirth 1999, Elbern et al. 1998] that it generally marks the 

14 boundary between tropospheric and stratospheric air. Poleward of the PV3.5 contours, 

15 there are relatively few clouds and FLO concentrations are much smaller. The PV3.5 

16 contours also enclose warm regions at 215 hPa, related to latent heat release from tropical 

17 convection, and low T values at 147 and 100 hPa, where adiabatic cooling in upwelling 

18 dominates as convective influence extends up to the cold tropopause. Higher 

19 stratospheric T values are found poleward of the PV3.5 contour. At 100 hPa, both MLS 

20 and GEOS-5 have an H 2 0 minimum over the Western Pacific, extending somewhat to the 

21 east of the lowest OLR and to the north of peak IWC, but coincident with minimum T. 

22 The convective regions over equatorial South America and Africa are warmer and 

23 moister at 100 hPa than that over the Western Pacific, consistent with the premise that T 
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1 controls humidity near the tropical tropopause [ Holton and Gettelman 2001; Read et al. 

2 2004], although convective dehydration may also play a role [Sherwood and Dessler 

3 2001], 

4 GEOS-5 IWC and PLO at 215 hPa are quite similar to MLS fields both in 

5 morphology and in magnitude, although GEOS-5 is moister than MLS at this level. 

6 GEOS-5 has less 147 hPa IWC and more H 2 O than MLS. This might be due to too much 

7 sublimation and/or too little condensation in the model's microphysics. GEOS-5 at 100 

8 hPa has smaller values of IWC and is drier in the tropics and wetter in the extra- tropics 

9 than MLS. The stronger latitudinal gradient of 100 hPa PLO in GEOS-5 may indicate 

10 some deficiencies in its representation of mass transport from the troposphere to the 

11 stratosphere. GEOS-5 is on average warmer than MLS in the tropics by ~3 K at 215 hPa 

12 and ~1 K at 147 hPa, which is not unexpected due to the known MLS low biases in 

13 cloudy regions (~2 K at 215 hPa and ~0.5 K at 147 hPa, Schwartz et al. [2008]). MLS 

14 and GEOS-5 tropical 100 hPa temperatures agree to within -0.5 K. 

15 Ligure 2 shows the 4-year (January 2005 to December 2008) tropical (15°S-15°N) 

16 mean profiles of IWC, H 2 O and T from both MLS and GEOS-5, along with their daily 

17 standard deviations from the 5-year mean. The GEOS-5 IWC profile agrees within 12% 

18 with MLS at 215 hPa to 147 hPa but becomes 30%, 50% and 70% smaller than the MLS 

19 IWC at 121 hPa, 100 hPa and 83 hPa, respectively. Although these are all within the 

20 estimated (factor of 2) uncertainty of MLS measurements, the smaller GEOS-5 IWC 

21 amounts above 147 hPa suggest that convection in the model does not extend to 

22 sufficiently high altitudes. 

23 The GEOS-5 and MLS differences in tropical H 2 0 and T appear to have a source 

24 other than that for the IWC differences. GEOS-5 H 2 O is limited to values corresponding 
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1 to 100% or less relative humidity, and an overestimate of T could possibly lead to an 

2 overestimate of HiO. Figure 2 shows that, after accounting for the known MLS ~2 K cold 

3 bias, GEOS-5 215 and 178 hPa tropical T is still larger than that of MLS by ~1 K, which 

4 may contribute to the large FLO in GEOS-5 at the two levels. However, the saturation 

5 HiO mixing ratio profile computed using GEOS-5 T suggests that this can not explain the 

6 HoO discrepancy between MLS and GEOS-5 at 147 hPa altitude and above. In summary, 

7 GEOS-5 215 hPa tropical H 2 O is larger than that of MLS by -50% and 215 hPa T is 

8 larger by -3.5 K, both of which are significant compared to estimated MLS measurement 

9 uncertainties. GEOS-5 and MLS tropical 147 hPa H 2 O agree within -30% (only slightly 

10 larger than the estimated 20% MLS measurement uncertainty) and 147 hPa T agree to 

11 -0.2 K (within the MLS measurement uncertainty). GEOS-5 and MLS tropical 100 hPa 

12 H 2 O agree to -15% and T to 0.2 K (both within the MLS measurement uncertainty). 

13 3.2. Seasonal maps 

14 Figure 3 shows seasonal MLS IWC and H 2 O maps at 100, 147 and 215 hPa. The 

_9 

15 overlaid 240-Wm ' contour generally encloses the highest values of both IWC at all three 

16 pressure levels and H 2 O at 147 and 215 hPa. High IWC in December to February (DJF) 

17 is concentrated south of the Equator in central-south Africa, the Western Pacific and 

18 South America. In June-August (JJA), the maximum IWC is distributed over the South 

19 Asian monsoon region, while South American convection has shifted northward to cover 

20 Central America. Seasonal variations over the Western Pacific are relatively small. The 

21 seasonal variation of the ITCZ (Inter-Tropical Convection Zone) and IPCZ (Inter-Pacific 

22 Convection Zone) is also apparent in MLS IWC. 

23 At 215 hPa, maxima in both IWC (Figure 3a) and H 2 0 (Figure 3b) are collocated 

24 with low OLR, indicating convective moistening of the UT in all seasons. At 147 hPa, the 
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1 HiO maxima are over the western Pacific in DJF, and over southern Asia in JJA, in both 

2 cases slightly north of the strongest convection. Studies using MLS data [e.g. Fu et al. 

3 2006; Park et al. 2007] have shown convectively-lofted FLO is trapped in the strong anti- 

4 cyclone over the Tibetan Plateau during the Asian summer monsoon, where IWC and 

5 FLO are seen distributed across the tropopause (PV3.5 contour) into the lower 

6 stratosphere. At 100 hPa, the minimum FLO values are found in the cold region over the 

7 tropical western Pacific in all four seasons. 

8 Comparing MLS and GEOS-5 maps (shown in Supplemental Figure 1), there is 

9 overall similarity in both IWC and FLO in terms of seasonal variations, with the 

10 differences shown in Figures 1 and 2 also evident in the seasonal maps. However, GEOS- 

11 5 shows a much less evident IWC ITCZ feature than MLS, especially in DJF and JJA. 

12 This is thought to be due to GEOS-5 underestimating the height of convective 

13 penetration, as has been mentioned earlier. We note, however, that GEOS-5 has less 

14 discrepancy with MLS than GEOS-4. 

15 4. Time evolution of UT H 2 O and IWC 

16 4.1. Latitude-time evolution 

17 Latitude-time sections of zonal-mean IWC and HoO from MLS (Figure 4a) and 

18 GEOS-5 (Figure 4b) further illustrate the seasonal evolution from August 2004 to July 

19 2009. The patterns of evolution of 215 and 147 hPa H 2 O and IWC, as well as 100 hPa 

20 IWC are qualitatively similar, while the 100 hPa H 2 O pattern is noticeably different in 

21 both MLS and GEOS-5. At 215 and 147 hPa in MLS and GEOS-5, the meridional 

22 movements of high IWC and HoO are in phase and follow the Sun, with highest IWC and 

23 H 2 O in the northern summer. While the GEOS-5 147 and 215 hPa IWC are in acceptable 

24 agreement (within -15%) with MLS, the GEOS-5 147 and 215 hPa H 2 O maxima are 
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1 larger than MLS by -30% to 50% throughout the year, as is also evident in Figure 1. 

2 MLS 100 hPa IWC is confined to a narrow latitudinal band, which shifts seasonally in a 

3 similar way to 215 and 147 hPa IWC. The underestimate of 100 hPa IWC in GEOS-5 

4 compared to MLS is clearly a year-round feature. The seasonal cycle in 100 hPa FLO is 

5 dramatically different from that at 147 and 215 hPa. The 100 hPa FLO is correlated with 

6 the annual cycle of temperature, with minima occurring over the Equator in boreal winter 

7 and spring. Annual maxima occur more or less simultaneously in both hemispheres 

8 (around August), with larger values of 100 hPa FLO in the northern hemisphere (NH) 

9 than in the southern hemisphere (SH). GEOS-5 NH (0°-60°N mean) 100 hPa H 2 O values 

10 are lower than MLS by -20%. In the SH (0°-60°S), GEOS-5 shows larger H 2 0 (-20% 

11 larger than MLS) which also persists longer throughout the year than in the NH. This is 

12 most likely due to problems in the model related to relaxation to a constant moisture 

13 mixing ratio in the stratosphere, which may not be related to the moist physical processes 

14 in the troposphere. 

15 There may be evidence for a two-year cycle in the magnitude of the H 2 0 minimum at 

16 100 hPa in both MLS and GEOS-5, which could be a manifestation of the tropical 

17 biennial oscillation (TBO) [e.g. Li et al. 2001]. Longer-term data are needed to determine 

18 if this is really the case. 

1 9 4.2. Height-time evolution 

20 Figure 5a shows the height-time section of tropical (15°S-15°N) daily mean MLS H 2 0 

21 anomalies from the 5-year (August 2004 to July 2009) mean, illustrating the so-called 

22 “tape-recorder” signal [Mote et al. 1996]. There is a clear vertical transport of H 2 0 from 

23 121 hPa through the stratosphere. An apparent two-year cycle in the magnitude of the 

24 tape-recorder signal, particularly evident in the magnitude of the boreal- winter minimum 
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1 around 100 hPa, may be related to the quasi biennial oscillation (QBO) [e.g. Baldwin et 

2 al., 2001] or TBO [e.g. Li et al. 2001]. H 2 O signals imprinted at the bottom of the 

3 stratosphere are maintained through the stratosphere for 12 to 18 months as the air rises. 

4 The tape recorder is less clear in the upper stratosphere, although the two intense dry 

5 phases are especially evident up to ~1 hPa. The seasonal cycle of tropical GEOS-5 H 2 O 

6 (Figure 5b) near 100 hPa has a similar magnitude to that of MLS. The amplitude of the 

7 GEOS-5 annual cycle at pressures below 68 hPa drops more rapidly than that of MLS, 

8 although the ascent rates are quite similar. This attenuation of the signal in GEOS-5 

9 arises because of its relaxation of stratospheric moisture to a constant value close to 6 

10 ppmv. Such a tape recorder signal does not appear in the IWC field, since the warmer 

11 stratosphere quickly sublimates ice particles, and ice is subject to sedimentation. 

12 4.3. Longitude-time evolution 

13 Figure 6 shows the longitude-time section of tropical (15°S-15°N) mean MLS IWC 

14 and H 2 O anomalies, relative to monthly mean data. On inter- annual time scales, El Nino- 

15 Southern Oscillation (ENSO) related signals dominate the variability. The El Nino (warm 

16 phase) patterns at 215 hPa (Figure 6b) are characterized by an enhancement of IWC and 

17 HiO in the eastern Pacific accompanied by a reduction of IWC and PEO in the western 

18 Pacific from late 2004 to early 2005 and late 2006 to early 2007. The opposite patterns 

19 are seen during the cold La Nina phase from late 2007 to early 2008 and late 2008 to 

20 early 2009. IWC anomalies at 100 hPa (Figure 6a) appear in phase with those at 215 hPa, 

21 but the PEO anomalies at 100 hPa are out of phase with the 215 hPa IWC and PEO. The 

22 PEO anomalies at 100 hPa are less localized in longitude than the IWC anomalies, and 

23 are strongest over the Indian Ocean (also see Figure 8 later in section 5). GEOS-5 IWC 

24 and H 2 0 data (see Supplemental Figure 2) show a similar pattern. 
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5. The UT response to ENSO 


2 The MLS simultaneous and collocated measurements of FLO. IWC and T provide an 

3 unprecedented characterization of the UT response to ENSO variations. Figure 7 shows 

4 time series of monthly mean Nino 3.4 SST (defined by Trenberth [1997] as the SST 

5 averaged for longitudes 170°-240° and latitudes 5°S-5°N), and the tropical (15°S to 

6 15°N) mean MLS IWC and FLO at 100, 147, and 215 hPa. The time series of tropical 

7 mean IWC and FLO anomalies have similar time evolution to that of the Nino 3.4 SST. 

8 We choose DJF 2005 and 2008 to represent, respectively, the warm and cold phases of 

9 ENSO. Figures 8a and b show the corresponding IWC, ILO and T anomalies for the two 

10 phases. A typical dipole pattern [ Semazzi and Indeje, 1999] is seen in IWC and FLO at 

11 215 and 147 hPa, with positive anomaly in the central-eastern Pacific and negative 

12 anomaly in the western Pacific during El Nino and the opposite pattern during La Nina. 

13 At 100 hPa in DJF 2005, a positive IWC anomaly and a negative PLO anomaly in the 

14 central eastern Pacific are accompanied by anomalies of opposite sign in the western 

15 Pacific. The opposite signatures are evident in DJF 2008. The negative IWC anomaly 

16 over the western Pacific during El Nino is an indication of reduced convection in 

17 response to warmer SST in the central-eastern Pacific [ Su et al., 2002], This is possibly 

18 associated with anticyclones west of the localized SST heating, as suggested by 

19 Highwood and Hoskins [1998]. The 100 hPa FLO anomalies are particularly strong in the 

20 Indian Ocean. Whether these particularly strong anomalies represent a teleconnection 

21 through “atmospheric bridge” [Alexander et al. , 2002; Klein et al., 1999] or a response to 

22 local SST anomaly is not clear. The ENSO response in GEOS-5 (see Supplemental 

23 Figure 3 and Figure 4) is similar to that of MFS. 
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1 6. Summary and Conclusions 

2 We have presented Aura MLS UT H 2 O and IWC measurements made from August 

3 2004 through July 2009, with comparisons to GEOS-5 analyses of these quantities for the 

4 same period. The global distributions of four-year-mean annual and seasonal averages, 

5 and tropical temporal evolution and response to ENSO, are given. Some comparisons 

6 with MLS and GEOS-5 UT temperatures are also discussed. 

7 Agreement between MLS and GEOS-5 H 2 O at 100 and 147 hPa is generally within 

8 the estimated MLS measurement accuracy of -20% (albeit slightly, but probably not 

9 significantly, worse, -30%, at 147 hPa). GEOS-5 has (during all seasons) smaller 

10 minimum tropical 100 hPa H 2 O values, and moister extratropics, than MLS, thought to be 

11 caused by the model’s relaxation to fixed stratospheric H 2 O concentrations. GEOS-5 215 

12 hPa H 2 O is larger than MLS by -50%, probably because the GEOS-5 deep convection 

13 does not extend sufficiently high. IWC agreement is within the factor of 2 estimated 

14 accuracy of MLS, but comparisons of IWC vertical distributions also suggest that GEOS- 

15 5 deep convection does not extend sufficiently high. There appears to be a significant 

16 difference in 215 hPa temperature, with GEOS-5 being -IK warmer after accounting for 

17 the known -2 K cold bias in MLS. MLS and GEOS-5 147 and 100 hPa temperatures 

18 agree on average to within -0.2 K, well within the MLS uncertainty. 

19 The tropical distributions of 215 hPa PLO and IWC are positively correlated; large 

20 values of both are associated with regions of deep convection, as previously found [e.g. 

21 Su et al. 2006a]. The distributions of 100 hPa H 2 O and IWC are negatively correlated, 

22 with less H 2 O and more IWC in regions of deep convection, as expected from “freeze- 

23 drying” of uplifted air. The transition from positive to negative correlation occurs 
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1 between 147 and 100 hPa. The tropical 215 hPa HoO and IWC seasonal variations track 

2 regions of deep convection, while the 100 hPa H 2 O seasonal variations track 100 hPa 

3 cold regions. The largest values of tropical HoO occur in the northern summer over the 

4 South Asia monsoon region; the smallest values of H 2 O occur in the northern winter over 

5 the western Pacific. 

6 Tropical zonal mean H 2 O and IWC exhibit strong seasonal and interannual variations. 

7 There is an indication of a two-year cycle in the tropical UT boreal winter H 2 O minimum, 

8 possibly related to the roughly two-year oscillation (TBO) in tropical SST [Li et al. 

9 2001]. The stratospheric H 2 O tape-recorder signal displays a two-year cycle, which may 

10 be related to the QBO [Baldwin et al., 2001]. GEOS-5 HoO appears to ascend faster 

11 through the upper tropical tropopause (121 hPa - 83 hPa,) than does MLS H 2 O, and has a 

12 smaller amplitude seasonal cycle in the stratosphere (where moisture is relaxed to a 

13 constant value). A future GEOS objective is to implement a more realistic stratospheric 

14 moisture module having methane oxidation chemistry. 

15 Fluctuations in tropical UT tUO and IWC are associated with moderate El Nino and 

16 La Nina events that occurred during the 5-year period analyzed here. PEO and IWC zonal 

17 mean fractional anomalies were -10%. The IWC and H 2 O deseasonalized 215 hPa 

18 anomalies exhibit a dipole pattern during El Nino (La Nina), with positive (negative) 

19 anomalies in the eastern Pacific and negative (positive) anomalies in the western Pacific. 

20 A strong positive (negative) 100 hPa H 2 O anomaly occurs over the Indian Ocean during 

21 El Nino (La Nina). 
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1 Figure Caption 

2 Figure 1: Annual mean IWC, PLO and temperature maps at 100, 147 and 215 hPa 

3 pressure levels from MLS observations (upper panels) and GEOS-5 analyses (lower 

_9 

4 panels). The black contour is the GEOS-5 OLR at 240 Wm . The grey contour is the 

5 GEOS-5 PV3.5. The GEOS-5 data are averaged onto 3°xl° boxes centered on MLS 

6 measurement locations. MLS averaging kernels are applied to GEOS-5 FLO and 

7 temperature data, and GEOS-5 IWC data are also vertically averaged in 3.5 km boxes 

8 centered on MLS data points. Four years of data from January 2005 to December 2008 

9 are used to compute the averages shown in this Figure. Thus this is a four-year average of 

10 “annual mean”. 

11 Figure 2: Tropical (15°S-15°N) mean IWC, H 2 O and T profiles from both MLS (black) 

12 and GEOS-5 (blue). The profiles are averages of daily mean tropical profiles from 

13 January 2005 to December 2008. The standard deviations of daily profiles for MLS 

14 (gray-shade) and GEOS-5 (blue-dashed) are also shown, as well as the saturation specific 

15 humidity profile (red) computed using the GEOS-5 Temperature. 

16 Figure 3: Seasonal mean MLS IWC (a) and FLO (b) maps at 100, 147 and 215 hPa 

_2 

17 pressure levels. The black contour is the GEOS-5 OLR at 240 Wm “. The grey contour is 

18 the GEOS-5 PV3.5. Data from December 2004 to October 2008 are used to compute the 

19 seasonal averages shown here. Each season includes 3-month from four different years. 

20 For example, JJA seasonal map is the average of June-August 2005, June-August 2006, 

21 June- August 2007, and June- August 2008; DJF seasonal map is the average of December 

22 2004-February 2005, December 2005-February 2006, December 2006-February 2007, 

23 and December 2007 -February 2008. 
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1 Figure 4: Latitude-time sections of zonal-mean IWC and H 2 O at 100, 147 and 215 hPa 

2 from (a) MLS observations and (b) GEOS-5 analyses, computed from daily zonal mean 

3 data. 

4 Figure 5: Height-time section of tropical (15°S-15°N) mean H 2 O anomalies from MLS 

5 (a) and GEOS-5 (b), computed from daily tropical mean data. MLS H 2 0 vertical 

6 averaging kernels are not applied to GEOS-5 H 2 0 in this plot. 

7 Figure 6: Longitude-time section of tropical (15°S-15°N) mean MLS IWC and H 2 O 

8 anomalies at 100 hPa (a) and 215 hPa (b), computed from monthly mean data. 

9 Figure 7: Time series of monthly mean SST (top panels) in the Nino 3.4 region 

10 (longitude 170°-240° and latitude 5°S-5°N), and (lower panels) monthly tropical (15°S- 

11 15°N) mean MLS IWC and H 2 0. 

12 Figure 8: Maps of 2005 DJF and 2008 DJF anomalies of IWC, H 2 O and T from MLS 

13 measurements. The anomalies are computed as the difference between 2005 or 2008 DJF 

14 averages and the four-year (2005-2008) average. Note that the 2005 DJF covers three 

15 months from December 2004 to February 2005, while the 2008 DJF refers to December 

16 2007 to February 2008 period. 
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Supplementary Figure 4 
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